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We present a macroscopic calculation of coherent electro-magnetic radiation from air showers 
initiated by ultra-high energy cosmic rays, based on currents obtained from Monte Carlo simulations 
of air showers in a realistic geo-magnetic field. We can clearly relate the time signal to the time 
dependence of the currents. We find that the the most important contribution to the pulse is 
related to the time variation of the currents. For showers forming a sufficiently large angle with the 
magnetic field, the contribution due to the currents induced by the geo-magnetic field is dominant, 
but neither the charge excess nor the dipole contribution can be neglected. We find a characteristic 
bipolar signal. In our calculations, we take into account a realistic index of refraction, whose 
importance depends on the impact parameter and the inclination. Also very important is the role 
of the positive ions. 



I. INTRODUCTION 

The main motivation for this work is the apparent 
need of radio detection experiments (LOPES 0,0], 
CODALEMA Q) for reliable calculations. Already in 
the earliest works on radio emission from air show- 
ers [1, 0, 01 ; a macroscopic treatment of the ra- 
dio emission was proposed, but at the time the as- 
sumptions about the currents were rather crude. In 
more recent work [1, 0, [l^, [ill, [l3|) microscopic 
picture of coherent synchrotron radiation from sec- 
ondary shower electrons and positrons gyrating in 
the Earth's magnetic field was employed. 

Recently, we performed macroscopic calculations, 
which allow one under simplifying conditions to ob- 
tain a simple analytic expression for the pulse shape, 
showing a clear relation between the pulse shape and 
the shower profile [l3|- The picture used was very 
similar to the one in Ref. [Bl, which has been refined 
by using a more realistic shower profile and where 
we calculate the time-dependence of the pulse. 

In the present paper, we advance further by com- 
puting first the four-current from a realistic Monte 
Carlo simulation (in the presence of a geo-magnetic 
field), and then solve the Maxwell equations to ob- 
tain the electric field, while considering a realistic 
(variable) index of refraction. Although this index 
varies only between 1 and 1.0003, this variation has 
quite interesting consequences. 

For the moment, we neglect the finite extension 
of the shower (pancake thickness and lateral exten- 
sion), at a given time. As shown in this is 
a good approximation for large impact parameters; 
for smaller ones the finite extension has to be consid- 
ered. There are no conceptual problems in extending 
the present approach to finite size sources, which is 
simply a three-dimensional integral over the point- 



like expressions. This will be discussed in a future 
publication. 

The shower moves with almost the vacuum ve- 
locity of light c. There is a constant creation of 
electrons and positrons at the shower front, with 
somewhat more electrons than positrons (electron 
excess). This is compensated by positive ions in the 
air, essentially at rest. The electrons and positrons 
of the shower scatter and lose energy, and there- 
fore they move slower than the shower front, falling 
behind, and finally drop out as "slow electrons / 
positrons". Close to the shower maximum, the 
charge excess of the "dropping out" particles is com- 
pensated by the positive ions, since there is no cur- 
rent before or behind the shower. Taking all to- 
gether we have the situation of a moving charge, 
moving with almost the vacuum velocity of Hght, 
even though the electrons and positrons are moving 
slower. For a consistent picture one must not forget 
the positive ions, which are usually not considered 
when talking about the shower. They are very slow, 
but their position of creation moves with the veloc- 
ity of light (shower front) , as well as the position of 
the "charge neutraHzation" at the "shower tail". The 
above picture is true with or without geo-magnetic 
field, in the former case the electrons and positrons 
move on curved trajectories due to the Lorentz force. 
Neglecting the finite dimension of the shower, one 
has a four-current 

j{t',x)^jit')s'{x-at')), (1) 

with a longitudinal component due to charge excess, 
and a transverse component due to the geo-magnetic 
field. The precise form of J{t') will be obtained from 
realistic shower simulations, using CONEX [l^. [l5|. 
Solving Maxwell's equations, we can express the 
electric field in terms of the four-current J and its 
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time derivative K = dJ/dt, evaluated at the re- 
tarded time. 

In this paper, we will develop expHcitly the for- 
maHsm of electric fields created from general point- 
like currents of the form shown in eq. |(T]), in par- 
ticular investigating carefully the role of a position 
dependent index of refraction. We analyze the elec- 
tric fields for some typical high energy showers, their 
dependence on the observer position and on the ori- 
entation of the showers with respect to the mag- 
netic field. We restrict ourself of sufficinently large 
impact parameters, so that singularities (related to 
Chrenkov radiation) do not matter. In all cases, 
we can clearly relate the electric field signal to the 
sources. The dominant contribution comes always 
from the time variation of the currents. For finite 
angles between the shower axis and the magnetic 
field, the biggest contribution is due to the time 
derivative of the transverse currents induced by the 
geo-magnetic field; but there is also a contribution 
(three times smaller, opposite polarization) due to 
the transverse current itself. There is also a sizable 
contribution from the time derivative of the charge 
excess, and from the dipole moment. The strength 
of these two contributions is roughly 1/4 each com- 
pared to the dominant contribution from the trans- 
verse current. The fact that the time variation of 
the currents is the dominant mechanism leads to a 
characteristic bipolar signal. 

These results actually confirm our earlier findings 
[l^ of a bipolar signal due to the time derivatives 
of the sources. This is in qualitative contrast to the 
findings of [1, [13, [ill, [l3|, where a microscopic 
picture has been employed, leading to a single pulse. 
The fundamental difference is that the time variation 
of the currents in the macroscopic picture treats in 
an effective way the creation and disappearance of 
charges, which contributes significantly. In a micro- 
scopic picture this has to be treated explicitly. 



II. GEOMETRY 

We characterize the trajectory of a moving point- 
like current (referred to as the shower trajectory) in 
the following way (see figE]): we imagine a sphere 
around the center M® of the Earth with radius 
i?0 -|- a, with i?0 being the radius of the Earth. The 
point A is by definition the intersection of the tra- 
jectory with the sphere (supposed the intersection 
exists). The angle between the trajectory and the 
axis M^A (referred to as vertical) is called 0. The 
trajectory is characterized by the point A, the an- 
gle (and an azimuthal angle, which is not impor- 
tant for the following discussion) . A position on the 
trajectory may be characterized by the distance d 




Figure 1: Characterizing the trajectory: A is the inter- 
section of the trajectory with a sphere of radius 7?® -I- a 
around the center of the Earth , with i?® being the 
radius of the Earth; the angle 9 measures the inclination. 



from the point A. Using these definitions, the height 
along the trajectory is given as 

h{d) = ^/{R(s + aY + d2 + 2(i?e + a)d cos 6 - i?©. 

(2) 

Let S be the point on the axis M®A representing 
the sea level. We define the "Earth frame" {S,Wx, 
Wy, Wz} (see fig. [21) such that Wz represents the 
vertical upward unit vector, and Wy is orthogonal to 
the "shower plane" defined by the shower trajectory 
and the point S. In other words, {S,Wx,Wz} is the 
shower plane. The shower velocity vector is given as 

(V sin 6 \ 
(3) 
—V cos 9 J 

with respect to the Earth frame basis. 




Figure 2: The trajectory as seen by an observer G situ- 
ated on the observation level. We use the Earth frame 
{S,Wx,Wy,Wz} to represent points and vectors. 

In addition to the sea level plane we have to con- 
sider an "observation level plane" at some height z 
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above sea level. We define an observation point G in 
the observation level plane with coordinates {x, y, z) 
with respect to the Earth frame. The point C rep- 
resents the center of the air shower front at time t\ 
given as SC = SA +Vt', with respect to the basis 
Wi. The vector R = CG is given as i? = SG — SC . 
We may define the point B of closest approach with 
respect to the observer. The distance from A is 



cLb — X said + (a — z) cos( 
The impact parameter is 



(4) 



h=\BG\^ \J [x cose- (a-z) sin0) +y2. (5) 
We also define for later use 



ts = ds/v 



(6) 



to be the time the shower passes at the position B. 

Due to the magnetic field of the Earth, electrons 
and positrons move on spiral trajectories with oppo- 
site orientations, representing a transverse current 
[1^ . This means, even when considering the shower 
as pointlike, the collective current J is not simply 
a product of charge times velocity, there is in fact 
a component transverse to the shower velocity. We 
compute this current by using the CONEX shower 
simulation program, as 



J 



E 



E 



(7) 



positrons i 



clccti' 



considering all electrons and positrons present at a 
given time. For calculating the currents, we employ 



6 'shower axis 




f ~ magnetic field 




at this point to have some idea about the times t' 
corresponding to shower maxima. The shower size 
and in particular the shower maximum depends pri- 
marily on the atmospheric depth, given as 



X = X{t') 



p{h{x)) dx, 



(8) 



integrating along the trajectory. Here, h is the 
height above sea level, and p is the air density. The 
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Figure 4: The depth versus the negative time —t', for 
different angles (using always a = 140 m). The horizon- 
tal line at X = 700 g/cm^ corresponds to the maximum 
of the charge distribution for showers with an energy of 
5 lO^'^eV. 

relation between X and t' (or d) depends of course 
strongly on the angle 6, as shown in fig. IH Since we 
are interested in negative time t' (before the shower 
reaches A), we plot X versus —t'. The horizontal fine 
corresponds to the maximum of the charge distribu- 
tion of showers with an energy of 5 lO^^eV, which can 
be used to obtain the times t' corresponding to this 
maximum. Varying the angle 9 from 0*^ to 90°, the 
times range from some 10 ps up to almost 2000 ^s. 
These numbers define the interesting time scales for 
the later discussions of retarded times. 



Figure 3: Reference frames. 

the shower frame {A,Ux,Uy,Uz}, see fig. [3l but the 
final results are expressed with respect to the Earth 
frame {SjWx, Wy, w^}. 

We will discuss the currents corresponding to air 
showers in more detail later, but it is useful already 



III. ELECTRIC FIELD DUE TO A 
GENERAL POINTLIKE SOURCE 

Let us consider a point-Hke four-current, 

j{t',x) = Jit')5\x-at')), (9) 

for a given trajectory ^{t'), corresponding to the 
shower position G in the notation of the preceding 
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chapters. A moving charge corresponds to = qc, 
J = qv, where c is the velocity of Hght. However, 
currents induced by the Earth's magnetic field re- 
quire a more general treatment, which we are going 
to present in the following. In this chapter we in- 
vestigate the case an index of refraction n equal to 
unity. 

We first define the basic quantities of our ap- 
proach. In order to obtain four-vectors x and 
we define a;° — ct and = ct' . We will need 
the time derivative of the four-current J and of 
the four-vector ^,and we define K = dJ/dx"^ and 
V = d£,{t') / dx'^ . We consider here the case of a con- 
stant velocity V . Finally, we define R — {x—^{t')) to 
be the four- vector joining the position of the current 
S,{t') and the position x of the observer. 

The four-potential at point x is found to be 



A 



Mo J 



47r \RV\ 



(see appendix) , where all quantities are meant to be 
evaluated at the retarded time t* . The three com- 
ponents of the electric field can be expressed as 



(11) 



Taking also the x dependence oft* into account, one 
gets (see appendix) 



1 



(12) 



47reoc \RV\^ ' 
with 

L'^f^ = RV R^kI^ - RV J'^V"" + VV R^jf^ (13) 

Concerning the denominator and possible singulari- 
ties due to its roots, it is useful to note that 



1 



\RV\ 



\dt*/dt\ 
\R\ 



(14) 



from which one can see that possible singularities of 
the field coincide with singularities of the derivative 
dt*/dt. 



IV. REALISTIC INDEX OF REFRACTION 

The index of refraction n in air varies spatially, 
mainly via its dependence on the altitude. We will 
study here the general case of an index of refraction 
being a function n{y) of the position y. In this case 
we define curvilinear coordinates via the light curves 
passing through the position y, see fig. [H The op- 
tical path length, s, is defined as the travel time of 
a light signal multiplied by c. An elementary path 




Figure 5: Light curves (full lines) in a medium with 
varying index of refraction. The dashed line represents 
a constant optical distance from the point y. 



length is expressed in terms of the components y^oiy 
as ds^ = n^dy^dy'' . The equation for the derivatives 
v'^ ~ dy'' /ds of the light curves is thus 

dv^ 



(10) ^ 



-n/i- 7 7 7 7 rt k 7 

= -T = — V V - 2 — V V 

ds ^ 



n 



n 



where we used 



n 



(15) 



(16) 



with rii = dn/dy^. The light curves can be ob- 
tained by integrating this differential equation start- 
ing from some position yb- One gets a family of 
curves, characterized by the initial velocity vectors 
vq (more precisely by the two angles characterizing 
the direction of the initial velocity, since the module 
is given). 

As the next step we discuss the effects of a varying 
index of refraction n{x) on the calculation of the 
electric field. The main difference compared to the 
case of n = 1 is the fact that the condition for the 
retarded time, expressed via 



(17) 



with = ct — ct' and R — x — £^, has to be modified 
by taking into account the curved light trajectories. 
We define a "distance" L = L{^,x) to he the optical 
path length from the source position £, to the ob- 
server position X, in other words L is the integral 
/ ds along the light curve from ^ to x. Instead of eq. 
(fl7| . we now have 



5{R°^-mt'),xf), 



which leads to 



A^ 



J 



with 



47r \RV\ 



R^^c(t~t*), R'^-L^L. 



(18) 



(19) 



(20) 
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All quantities are meant to be evaluated at ^(t*), 
with the retarded time t* corresponding to the roots 
of the argument of the 6 function in eq. ifTSl) . We 
thus get formally the same results as in the case 
n = 1, but with R replaced by R, and a different 
definition of t* . The computation of the retarded 
times will be discussed in the next chapter. 

Taking the derivative of i?° — L'^, we get the 
useful relation 



R" 
RV' 



with 



= c(t - t*), R' = L—L 



(21) 



(22) 



which again relates possible singularities of the po- 
tential A (roots of the denominator) to singularities 
ofdt*/dt. 

To compute the derivative d^A^ of the potential, 
we use 

d'^RV ^V" -VV d^'ct* (23) 
with the four-vectors V and V given as 

v'^^r^'v^, v-'^r^Vf,, (24) 

where the tensors g and g are defined as 

(25) 



d d d d 

Making use of eq. (|23l) together with eq. l(2T|) . we 
finally can write the electric field as 



1 L 



iO rOi 



47reoc \RV\3 ' 



(28) 



with 



L"'^ = RV R"K^ - RVV°'J'^ + VV R^jf^. (29) 

In case of n = 1, we have R = R = R, as well as 
9 — 9 = 9t cind therefore V — V — V , and we recover 
the formula derived earlier for the case n = 1. 

It is of course not surprising to rediscover the 
n — 1 case, but it is important to notice that the 
equation for the the field E'^ in case of n > 1 has 
the same structure as the one for n = 1, just with 
" ~ " and " ~ " decorated quantities instead of the 
"bare" ones. The former ones depend continuously 



and smoothly on a parameter i5 which characterizes 
the deviation of n from unity. We have, for exam- 
ple, R = R{S), with Hm^^o ^ — R, and thereafter 
R = R + i?(i)(5 + .... Since 5 < 0.0003 we can ac- 
tually already ignore to a good precision the linear 
term, in other words we can forget about all the "~" 
and "~" symbols, with two exceptions: first the term 
\RV\ in the denominator can become very small, and 
here even a very small value of S could matter. Sec- 
ond the term VV is certainly very small {VV — 0), 
but can be comparable to RV, close to a singularity. 
We finally obtain 



E' 



1 



47reoC \RV\^ 



(30) 



with 



L""!^ = RV R^Rl^ - RVV^J^^ + VVR''J'^. (31) 
In vector notation, we may write 



E — Erk + Ejv + Erj 



with 



_ Wrk p 



w. 



JV 



E 



D ' D 
using D = AireQclRV}'^ , and with 

Wab = AB° - A°B 



RJ 



(32) 



VV Wrj 

RV D ' 
(33) 



(34) 



Here, we assume RV > 0, to be verified later. The 
first two terms in eq. ( [33|) are the contribution Ejv 
due to the current J, and the contribution Erk due 
to the variation K — c^^ dJ/dt' of the current. The 
third contribution Erj is small, unless one is close 
to a singularity, due to the smallness of the ratio 
VV/RV, for high energy showers with \V\ « 1. 

As already mentioned (eq. l(2T|) ). \RV\ may be 
expressed in terms of the derivatives of the retarded 
time, as \RV\ = c(t - t*) / {dt* / dt). We there- 
for should be very careful about the evaluation of 
the retarded time t*{t) for a given observer time t, 
where the precise values of n may be important even 
though n — 1 is very small. 

V. RETARDED TIME 

The retarded time is obtained from solving L = 
i?", which may be written as 



L{i{t'),x) - c{tB - t') ^ c{t ~ te), 



(35) 



where a term —cts has been added on both sides, 
for later convenience. The time Ib corresponds to 
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X 



Figure 6: Light curve joining x and ^. 

the time when the shower is closest to the observer 
position X. 

We first evaluate the l.h.s. of eq. (|35|) . for a given 
arbitrary value of t' , corresponding to some shower 
position £,{t'). To obtain L, we solve the differential 
equation 

^ = !^,V-2%V, (36) 
as as n n 

with the boundary conditions 

y(0) = f, yis*)^l (37) 

with s* to be determined (see fig. [6]). The index of 
refraction n is 

cm"^ 

n(y) = 1 + 0.226 p{h{y)), (38) 

9 

with p being the density of air and h{y) the at- 
mospheric height h corresponding to a position y. 
We employ an adaptive step-size BuHrsch-Stoer al- 
gorithm to integrate very efficiently the differen- 
tial equation, and the three-dimensional Newton- 
Raphson formula to accommodate the boundary 
conditions. Having solved this problem means hav- 
ing found an s* which corresponds to the optical 
path length joining ^ and x, so 

mt'),x) = s*. (39) 

To understand the results, it is useful to consider 
two limiting cases, namely n — 1 and n = riground ~ 
1.0003 (the reality lying in between). For constant 
n, the l.h.s. of eq. l(35l) reads 

L-c{tB-t')=n\R\-c{tB-t'), (40) 

where R has been defined earlier to be the vector 
joining the emission point on the shower and the 
position of the observer. One gets 

L - c{tB - t') = ny/v^{tB -t'Y + b'^ - c{tB - t'). 

(41) 

The expression L — c{tB — t') divided by b is an 
universal function oi t = c{tB — t')/b, namely 

nV/3V2 + 1 - r, (42) 



with l3 = v/c. For n = 1 and /3 < 1, we have 
a monotonically decreasing function with no lower 
limit. In all the following considerations, we always 
consider a fast shower, with v = c, so f3 = 1. In 
this case, we still have a monotonically decreasing 
function, but with a lower limit equal to zero (for 
T oo). Keeping (3=1, but taking n > 1, the curve 
reaches a minimum at some finite value of t, and 
then increases. We show the corresponding results 
in fig. [71 where we plot L — c{tB — t') versus ts — t' , 
for two impact parameters b and several inclinations. 
The dashed curve refers to n = 1, the dotted one 
to n > 1. To underline the scaling behavior, we 
scaled abscissas and ordinates as 1/&, and indeed we 
observe, for all angles and impact parameters, the 
same curves, for the two limiting cases with constant 
n. 

Also shown in fig. [7] are the curves L — c{tB — t') 
versus ts — t' , for the case of a realistic index of 
refraction (full lines), obtained from the numeri- 
cal procedures discussed above. As expected, these 
curves are well in between the two limiting cases 
n = 1 and n = nground- The scaling is violated: de- 
pending on impact parameter and inclination, the 
exact curve follows more closely either the n = 1 or 
the n = Uground case. Whereas for vertical showers 
the exact curves are close to the n = 1 case, we see 
an increasing deviation with increasing inclination, 
up to the extreme at 90°, where the curves are close 
to the results for n = riground- All this is understand- 
able, since, with increasing angle, even light emitted 
at early times (large ts —t') travels more and more 
at n « nground- It should be noted that at very large 
tB—t', the curves L — c{tB—t') approach a constant , 
for obvious reasons, and there will be a (at least one) 
turning point. 

The retarded time t' = t* is the solution of 

L - c{tB - t') = c{t - tB), (43) 

for a given observer time t. In other words, t* is the 
intersection of the curves L — c{tB—t'), shown in fig. 
Ill with a horizontal line with ordinate c{t — Ib)- In 
the range oUb — t' up to 3000 /iS (useless to go to 
earlier times), we have no solution for t < <min, while 
for larger values there are one, two, or even more so- 
lutions. For the minimum time we have tmin > ts, 
where the precise value depends on angle and im- 
pact parameter. We compute the retarded times t* 
numerically. A first estimate is obtained via inter- 
polation using tables of L — c(is —t'), a refinement 
is obtained by again solving the differential equation 
eq. (|36| . but now with the boundary condition 

m=S, y{t~t*)=C (44) 

Only few Newton-Raphson iterations are needed to 
get a very high precision. The results are shown in 
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1 10 10^ 10^ 1 10 10^ 

tg-t' {[is) tg-f {[is) 



Figure 7: The optical path length L{^(t'),x) minus c(tB — t') versus ts — t' , for impact parameters of 500m and 
100 m, for different angles . The solid lines represent the results for a realistic variable index of refraction. The other 
curves refer to calculations for a constant n: namely n = 1 (dashed) and n — rigiound ~ 1.0003 (dotted). 



fig- El where we plot ts — t* versus t — ts, with t* 
being the retarded time. We get of course what we 
expect after the above discussion. In case of n = 1, 
we have no solution for negative t — ts, for positive 
values there is exactly one solution, representing a 
monotonically decreasing curve. The other limiting 
case n = Aground provides a two- valued curve and we 
have two solutions for t — ts being bigger than some 
minimum value, no solutions otherwise. The results 
for a realistic index of refraction (full lines) lies in 
between the two limits. For large impact parameters 
(like 500 m) , the exact curves are close to the n — 1 
case, with increasing angle and decreasing impact 



parameter, the curves get more and more close to 
the other limiting case, leading to two or even more 
retarded times, for a given observer time t. This 
means that for the electric fields one has to sum up 
the contributions due to several retarded times, if 
the currents J{t*) at these times are non-zero. 

From the above discussion - and in particular from 
the comparison of the exact calculation with the lim- 
iting cases of constant index of refraction - one gets 
the impression that the main effect is not really due 
to curved light trajectories but to the fact that the 
index changes along the trajectory. In order to verify 
the conjecture, we compute the optical path length 
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Figure 8: The retarded time expressed via —t* versus observer time t — ts, for impact parameters of 500 m and 
100 m, for different angles . The solid lines represent the results for a realistic variable index of refraction. The other 
curves refer to calculations for a constant n: namely n = 1 (dashed) and n = riground ~ 1.0003 (dotted). 



in "straight line approximation" as 

Lstrmghti^{t'),x) ^ / n{y)ds, (45) 
Jgc 

where GC is the straight Une joining abserver posi- 
tion X and charge £,{t'). Subsequently we compute 
the retarded times by solving Lstiaight — R^- The 
results are shown in fig. [9] as dotted lines, together 
with the exact curves (full lines) . The two results are 
almost indistinguishable, with the only exception of 
a small deviation for b — 500 m and — 90°, at very 
early retarded times. The straight line approxima- 
tion is therefore for most applications largely suffi- 



cient. 

A delicate point is related to the derivative dt*/dt. 
We have 



RV C{t-t*) dt' ^ 

showing that possible singularities of the expression 
for the electric field, corresponding to the roots of 
RV , are related to singularities of dt* /dt. In the 
simplest case of a large impact parameter and small 
angles, is — i* is a monotonic function of t — ts, 
and thus dt* /dt is as well, as shown in fig. [Tol 
Compared to the case n — \ (dashed) the exact 
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Figure 9: The retarded time expressed via ts —t* versus observer time t — ts, for impact parameters of 500 m and 
100 m, for different angles, considering a realistic variable index of refraction. The solid lines represent the exact 
results, the dotted lines refer to the "straight line approximation". 



curve (full line) grows more rapidly for t — ts ^ 0- 
For larger angles, we obtain singular points, where 
dt* /dt diverges. This means that whenever we have 
a situation where dt* /dt is singular, and in addition 
the current is non- vanishing at the corresponding re- 
tarded time, we have to consider the finite extension 
of the source, which amounts to an integration over 
the (square root) singularity, which gives a finite re- 
sult. 

It should be noted that the singularities of dt* /dt 
correspond to the singularities that gives rise to 
Cherenkov radiation. It refiects the fact that elec- 
tromagnetic waves emitted at different times arrive 



simultaneously at the observer position. 

The retarded time where the first singularity may 
occur depends on the impact parameter: it is pro- 
portional to b. It is somewhat larger than 10 /xs for 
b — 100 m, corresponding to a source more that 3 km 
away. 

To avoid this discussion related to singularities for 
the moment (to be treated in a forthcoming publi- 
cation) , we restrict ourselves to situations which are 
"well behaved" in the sense that dt* /dt is nonsingu- 
lar there where the currents are non- vanishing. 
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electric field is given as a sum of three terms, 



5 

* 

-a 




t - 1 ()is) 



Figure 10: The derivative dt*/dt versus t — Ib, for 
h — 500 m and 9 — 27", in case of a realistic index of 
refraction (full line) and the case n — 1 (dashed). 
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Erj — 



VV_ Wrj 
RV 



D ' 
(47) 

with Wab = AB" - A^B. We may split the vectors 
W into components parallel and orthogonal to the 
shower axis. The components of the shower velocity 
are — 1, V-^ = 0, V^^ = (3, with (3 being very 
close to unity. We thus find 



W 



JV 



J°(l-/3) wO 



and 



(48) 



(49) 



So the field Ejv is almost purely transverse, parallel 
to the transverse current, as sketched in figs. [TT] 
, where we show the electric field vectors for four 
different observers situated in a plane transverse to 
the shower axis. 



VI. THE ROLE OF THE POSITIVE IONS 



N 



For the individual charged paticles the magnitude 
of turns out to be always bigger than that of 
J^, even if one turns off the magnetic field. This 
is due to the fact that the electrons move slower 
than the shower front. They are constantly loos- 
ing energy and eventually drop out. Close to the 
shower maximum we have a steady state since here 
the number of net charge created is equal to the net 
charge dropped. On the other hand, the net charge 
created must correspond to a corresponding posi- 
tive charge in terms of ionized atoms, and therefore 
the net charge behind the shower (and in front of 
it anyway) is zero. So we have effectively a moving 
charge, moving with the speed of the shower front, 
which means that the z and the component of the 
total, collective, effective current must be equal. As 
mentioned above, we do not get this identity if we 
would only consider electrons and positrons, and not 
the ions. The ions will contribute to J", not to the 
other components of the current (they are not mov- 
ing) to yield finally J^^^ = J^^^ = J^, where "tot" 
refers to the fact that we consider all, electrons and 
ions. In the following calculations, we will always 
use J°Qt — , although strictly speaking this is only 
vaHd close to the maximum, but on the other hand 
the most important contribution to the fields is due 
to the charges close to the maximum. For the rest of 
the paper J will refer to the total current, including 
the ion part. 

The property J° = has important conse- 
quences concerning the polarity of the signal. The 



W transv. E 
current 



Figure 11: The electric field vectors from the contribu- 
tion Ejv, being proportional to the transverse current 



Furthermore, we find W^j^ « RV K° , (where 
here means correct to a precision of n — 1) which 
can be proven by using 1 = {d^ — 9ll)ci*« (i?" — 
R^/RV, due to time translation invariance. Since 
the field in general is large when RV is small, we 
may conclude that wj|^ is small compared to the 



transverse component, which we express as 



wIk = K° 



For the transverse part we find 



RK 



R^K'^-R'^K^, 



(50) 
(51) 



which means that here not only the transverse cur- 
rent (or better: its time derivative) contributes, but 
also the longitudinal one via Whereas the lat- 
ter one is radial with respect to the axis, the former 
one is anti-parallel to the direction of the transverse 
current, see fig. [HI 

We said earher that the third term Erj is small, 
away from singularities; let us investigate this term 
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Figure 12: The electric field vectors from the contri- 
bution Erk; one part is proportional to the time varia- 
tion of the transverse current J^(left plot), a second 
part is proportional to the variation if " of the longitudi- 
nal current, and proportional to i?^ and therefore radial 
(right plot). 



more closely. Using the same arguments as for Wrk , 
we find 



VV RV J° 



7° 

yyl. 



0, (52) 



where 0" means small compared to the transverse 
component, which is 



E 



RJ 



RV 



D 



(53) 



Lets us consider (just here, for an estimate) a con- 
stant index of refraction n > l^and^let us assume 
that p > n.-^. The prefactor VV/RV would then 
be proportional to 7~^(1 — nf3 cosa)~^ , with a be- 
ing the angle between the shower axis and R. Even 
though 7 is large, the term contributes when n(3 cosa 
approaches unity, which corresponds to the above- 
mentioned singularity. This is actually the well 
known "Cherenkov radiation". It should be noted 
that these singularities have to be treated very care- 
fully, using the realistic position dependent index of 
fraction n{y), and naively replacing V RhyV and R 
is certainly not permitted. In this paper we restrict 
ourselves to not too small impact parameters, where 
the contribution Erj is negligible. 



VII. DIPOLE FIELD 

So far we discussed electric fields created by the 
current due to the charge excess - taking into ac- 
count the positive ions - and the fields due to 
transverse current induced by the geomagnetic field. 
However, this magnetic field does not only create a 
transverse current, but also a dipole moment P orig- 
ina ting from the separation of the electric charges 
[1, We define q to be the product Ne of the 

number of electrons and positrons multiplied by the 



elementary charge e, and we define the derivatives 
of q as qn = c~"'d"'q/ dt'" , for n = 1,2. The ratios 
— jq are slowly varying (compared to the 
variation of q) along the trajectory. The potential 
due to this dipole moment may be written as 



(54) 



where A^^ is the potential due to the current Jq — 
cqV, representing a moving charge q. The electric 
dipole field is then given as 



E. 



dip 



'M'^dkE'^, 



(55) 



with Eg being the electric field due to a moving 
charge q. We may use the formulas derived earlier, 
where the electric field is given as a sum of three 
contributions. Only one term (Erk) is dominant, if 
we are not close to a singularity, with K = cqiV. 
The dipole field can then be written as 



M'^dk qi 



R'V° ~ R^V^ 



M^q2{dkct*) 



+qi dk 



(i?V)2 

R'V° - R°V' 
(RV02 



The identities 



dkRV = 

dkct* = 



= g,^-V^dkct* 

Vk^VVdkct*, 

Rk/RV, 



(56) 
(57) 

(58) 



(59) 
(60) 
(61) 



allow us to obtain 

ill 

-917^+92 



{MR) Wrv 2 VV (MR) Wrv 
■+qi — 



with the denominators given as 



47reo (W)", 



(62) 



(63) 



and withW^B = AB° - A°B. The last term in eq. 
((62l) only contributes close to a singularity, and for 
all the applications in the following chapters we only 
consider the first two terms. 

We may again split the vectors into longitudinal 
and transverse parts, with respect to the shower 
axis. Let us consider the dominant second term. We 
have Wrv ~ 0, following the same arguments which 
lead to the conclusion that Wrk is small. Therefore 
also the dipole field is essentially transverse. Using 
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Figure 13: The electric field vectors from the dominant 
dipole contribution, being proportional to the projection 
of M to the radial vector pointing to the observer. 



W^y — i? and the fact that M is purely transverse, 
we get 



92- 



(MR^) R^ 



(64) 



which means that the field is proportional to the 
projection of M to the radial vector pointing to the 
observer, as shown in fig. [I3l 



VIII. SHOWER PARALLEL TO THE 
MAGNETIC FIELD 

In the following, we consider the case of a shower 
with an initial energy of 5 • 10^'' eV, an inchnation 
9 = 27", and an azimuth angle ip = 180°, defined 
with respect to the magnetic north pole, see fig. [141 
The angle ip refers to the origin of the shower, it is 
related to the angle (p from the preceding chapter as 
tp = 180° - ip. So i; = 180° thus imphes that the 
shower moves from south to north. We consider the 
magnetic field at the CODALEMA site, i.e. \B\ = 
47.3/xr and a = 153°, so the shower moves parahel 
to the magnetic field. We consider an observer at 
X = 0, y = 500 m, z — 140 m. We suppose a — z 
(so the shower hits the ground at 2: = y = 0). The 
coordinates x, y, z are defined in the Earth system. 




Figure 14: South- north shower (61 = 27°, ip 
being parallel to the geo-magnetic field. 
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Figure 15: The number A'^ of electrons and positrons, as 
a function of the time ts — t ■ 
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Figure 16: The x- (full line), y- (dashed) and z- 
component (dotted) of of the current J divided by Nec 
(upper plot), and the x- (full line), and y-component 
(dashed) of of the dipole moment P divided by A'^e, 
(lower plot), as a function of the retarded time. The 
vector components refer to the shower frame. 



The electron number N reaches its maximum at 
about — 15/iS, reaching almost 4 • 10^ particles, as 



shown in fig. [TH It is useful to plot the 

components of the current, divided by Nec, with 
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Figure 17: The retarded time expressed via ts— t* versus 
observer time t - ts, for b = 500m and 9 ^ 27° . The 
solid lines represent the results for a realistic variable 
index of refraction, the two other curves to n = 1 and 

^ Aground • 



e being the electron charge, and c the velocity of 
light, see fig. [TH The components J^/Nec and 
jy /Nee represent the average transverse drift ve- 
locity, in units of c, caused by the magnetic field. 
For the present example, the components /Nec 
and jy /Nec are zero. This is understandable, since 
the average transverse drift velocity is zero for show- 
ers moving parallel to the magnetic field. The com- 
ponent /Nec represents the charge excess which 
stays almost constant at a value around -0.2. The 
time variation of the currents is mainly due to the 
time variation of N . For the numerical calculations, 
also the small time variations of /Nec are taken 
into account. In fig. [TBI we also plot the components 
and py of the dipole moment divided by A'^e. As 
expected, the dipole moment is close to zero. It is 
not exactly zero due to statistical fiuctuations in re- 
gions where Ne is small. The actual Monte Carlo re- 
sults are strongly fiuctuating around zero, the curves 
shown in fig. [16] is obtained after smoothening the 
"raw data". 

In the expressions for the electric fields at some ob- 
server time t, the currents J and their derivatives K 
appear, evaluated at the retarded time t* . Following 
the procedure discussed earlier, we compute t*[t), as 
shown in fig. [iTl We see that a shower maximum 
at Ib — t* ~ 15 /iS corresponds to t — ts ~ 0.1 /iS. 
Given the time dependence J{t') of the current and 
the relation t* [t) between the retarded time and the 
observation time t, we may present the current as 
J{t*{t)) versus the observer time. In fig. [181 we show 
the components of the current and their derivatives 
in the Earth frame. The y— component vanishes, as 
expected, since this is the component perpendicular 
to the shower plane (containing the shower axis and 
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Figure 18: The x- (full line), y- (dashed), z- (dotted), 
and 0-component (dashed-dotted) of the current J and 
their derivatives (bottom) versus the observer time t. 
Here we consider the total current, including the ion 
contribution. The vector components refer to the Earth 
frame. 



the vertical axis). One should note that the width 
of J{t') is of the order of 10 /iS, whereas the width of 
J versus observer time is just a fraction of a micro 
second. 

We are now ready to discuss the electric fields. As 
shown earlier, we have 



E — -Edip + EjiK + EjVi 
with the dipole field i?dipand with 

Wrk ^ Wjv 



E 



RK 



D 



E 



JV 



D 



(65) 



(66) 



usmg 



and 



D = ATreon^iRV)^ 



Wab = AB° - A^B, 



(67) 
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Figure 19: The x- (full line), y- (dashed) and z-component (dotted) of Wrk (a), Wjv (c), and of the corresponding 
fields Erk (b), Ejv{d), the dipole contribution Edip (e), and the components of the total electric field, with all 
contributions (f). The vector components refer to the Earth frame. 



for two four- vectors A and B. The three compo- 
nents of the vectors W and the corresponding field 
contributions, as well as the three components of the 
complete field E are shown in fig. [191 The dipole 
field is of course negligible. The contribution Erk - 
having its origin in the time variation K of the cur- 
rents - is by far dominant. Only the y— component 
is sizable, meaning that the field is oriented radially 
compared to the shower axis (the observer sits on 
the y— axis). The contribution Ejy is very small, 
in comparison. This can easily be understood since 
here, without transverse drift, we have essentially 
a longitudinally moving charge, and in such a case 
JV" -J°V vanishes. 



IX. SHOWER WITH A LARGE ANGLE 
WITH RESPECT TO THE MAGNETIC 
FIELD 

In the following, we consider a shower with an 
initial energy of 5 • lO^^eV, an incHnation 6 = 27°, 
and an azimuth angle = 0°, defined with respect 
to the magnetic north pole, see fig. [201 

The angle ip refers to the origin of the shower, it 
is related to the angle (p from the preceding chapter 
as ip = 180° — ip. So ■0 = means the shower moves 
from north to south. We consider again the mag- 
netic field at the CODALEMA site, so the shower 
makes an angle of roughly 54" with the magnetic 
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Figure 20: North-south shower {9 = 27° , ip = 0°^ 
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Figure 21: The x- (full line), y- (dashed) and z- 
component (dotted) of of the current J divided by Nec 
(upper plot), and the x- (full line), and y-component 
(dashed) of of the dipole moment P divided by A'^e (lower 
plot), as a function of of the time Ib — t' . The vector 
components refer to the shower frame. 



field. We consider an observer at x — 0, y — 500 m, 
z = 140m (which means: east of the impact pent). 
We suppose a = z (so the shower hits the ground at 
X = y = 0). We plot again the components of the 
current, divided by Nec, with e being the electron 
charge, and c the velocity of light, see fig. [2TJ The 
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Figure 22: The x- (full line), y- (dashed), z- (dotted), 
and 0-component (dashed-dotted) of the current J (top) 
and their derivatives (bottom) versus the observer time 
t — ts- Here we consider the total current, including the 
ion contribution. The vector components refer to the 
Earth frame. 



components /Nec and /Nec represent the av- 
erage transverse drift velocity, in units of c, caused 
by the magnetic field. Contrary to the case stud- 
ied in the preceding chapter, we have now a non- 
vanishing transverse drift and /Nec is non-zero. 
The other component /Nec vanishes because both 
the shower trajectory and the magnetic field are in 
the y — z plane (in the shower frame). The compo- 
nent /Nec represents the charge excess. The di- 
mensionless quantities /Nec vary only little with 
time, corresponding to a transverse drift velocity of 
0.04 and a fractional charge excess of -0.2. The small 
difference compared to the result from the previous 
chapter is actually due to event by event fiuctua- 
tions. The time variation of the currents is mainly 
due to the time variation of N . For the numerical 
calculations, also the small time variations of /Nec 
are taken into account. Not only the transverse cur- 
rent is non- vanishing, also the dipole moment con- 
tributes, see fig. [2T1 

Expressing the retarded time in terms of the ob- 
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Figure 23: The x- (full line), y- (dashed) and z-component (dotted) of Wrk (a), Wjv (c), and of the corresponding 
fields Erk (b), Ejv{d), the dipole contribution Edip (e), and the components of the total electric field, with all 
contributions (f). The vector components refer to the Earth frame. 



servation time, the time dependence J{t') of the cur- 
rent can be expressed asj(t*{t)) versus the observer 
time. In fig. [221 show the components of the 
current and their derivatives in the Earth frame. 
We now discuss the electric fields, 

E = Edip+ERK+Ejv - ^dip + -^ + -^. (68) 

The three components of the vectors W and the cor- 
responding field contributions, the dipole field com- 
ponents, as well as the three components of the com- 
plete field E are shown in fig. [231 The largest contri- 
bution is Erk , which has its origin in the time vari- 
ation of the currents. Also Ejv also contributes sig- 



nificantly. This contribution is broader but smaller 
than Erk, and of opposite sign, so the sum of the 
two gives finally a bipolar field. Compared to the 
shower parallel to the magnetic field, studied in the 
previous chapter, both contributions are much big- 
ger, which means that the dominant contribution 
is due to the transverse current caused by the geo- 
magnetic field. And here, the contribution due to 
the time variation of this current is by far domi- 
nant. But we learn as well that the charge excess 
contribution cannot be neglected. Also the dipole 
field contributes significantly, with an opposite sign 
relative to Erk- 
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X. DIFFERENT OBSERVER POSITIONS 

In this chapter, we will compare the fields seen by 
an observer at four difi^erent positions: south, east, 
north, and west of the impact point, always the same 
distance (500 m). We show in fig. [HI the retarded 
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Figure 24: The retarded time expressed via fs— t* versus 
observer time t — ts, for a shower with 9 — 27° , tp — O'^ , 
seen by observers south (full ), or east (dashed) of the 
point of impact, at a distance of 500 m. 

time versus observer time t — ts, for our shower with 
9 — 27° , i/" = O", as seen by different observers. 
For a given 6, the curves are different for different 
observers. This fact will help to understand the dif- 
ferent widths of the signals seen from the different 
observers. 

We consider first of all a vertical shower {9 = 0°), 
with the magnetic field artificially turned off. South 
corresponds to positive x, east to positive y. The 
results are shown in fig. [251 We observe a radial 
electric field, as expected. The sHght deviation of 
a perfect axial symmetry is due to the fact that the 
calculation is based on a Monte Carlo generated cur- 
rent, which is of course only approximately symmet- 
ric. So even without geo-magnetic field, one observes 
a sizable field of the order of 25 /xV/m. 

Also shown in fig. [25] are the fields for a slightly 
inclined shower, with 9 — 27°, %p — 0°, still with- 
out geomagnetic field. The results are quite similar, 
we still observe a radial electric field, with respect to 
the shower axis (with gives a non-zero z— component 
in the Earth frame in case of a north or south ob- 
server). The signal is also narrower, which is also 
easy to understand, since in case of non-zero incli- 
nation the maximum of the current corresponds to 



an earlier retarded time t* compared to a vertical 
shower, and correspondingly the derivative dt* /dt is 
bigger. Therefore a current with a given width in t* 
will be narrower in t, in case of an inclined shower, 
leading to a narrower signal. 

In fig[26| (upper four plots), we show the results 
for the same incHned shower, with a reaHstic geo- 
magnetic field, but without the dipole contribution. 
We get in addition to the charge excess contribution 
(seen in the case without magnetic field) a strong 
additional contribution perpendicular to the plane 
defined by shower axis and geo-magnetic field. As 
discussed earher, the dominant contribution is E^k^ 
due to the time variation K of the currents, and here 
mainly due to the transverse current from the geo- 
magnetic field. However, seen from the comparison 
with the B — Q case, the charge excess is not negli- 
gible. 

We finally show in fig. [261 (lower four plots) the 
complete fields for the four observers, including the 
dipole field. Whereas for the north and the south 
observers including the dipole field does not change 
the results, we find for the east and west observers 
a considerable reduction of the strength due to the 
dipole contribution. 



XI. CONCLUSION 

We introduced a macroscopic calculation of radio 
emission from air showers initiated by ultra-high en- 
ergy cosmic rays, based on currents obtained from 
CONEX Monte Carlo simulations of air showers in 
a reaHstic geo-magnetic field. We took into account 
a realistic index of refraction n, which varies along 
the trajectory of the shower. 

We demonstrated that in the case of a realistic 
index of refraction n{y) (depending on the space 
position y) one can derive expressions for the elec- 
tric field which are formally equivalent to those for 
n = 1, however with a modified definition of the 
retarded time. One has to take into account the 
fact that light moves in general on curved trajec- 
tories, and that the velocity of light changes along 
the trajectory. These curves are solutions of differ- 
ential equations containing n{y). Finding the re- 
tarded time amounts to integrating this differential 
equation, with given boundary conditions. An ac- 
curate approximation can be obtained by assuming 
straight fight trajectories. Another deHcate point 
is denominator D in the expressions for the electric 
field, which is inversely proportional to the square of 
the time derivative dt* /dt of the retarded time. This 
derivative may become singular, and the position of 
the singularities is strongly affected by the precise 
value o{ n. So even though n is very close to unity. 
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Figure 25: The x- (full line), y- (dashed) and ^-component (dotted) of the electric field of a vertical shower (upper 
four plots) and a shower with 6 — 27", ip = O" (lower four plots), without geo-magnetic field, for four observers 500 
m south, east, north, and west of the point of impact. The vector components refer to the Earth frame. 



the small deviations from unity are important. 

We showed that the electric field can be expressed 
in terms of the four-vector R joining the observer's 
position and the shower, the velocity four-vector V 
of the shower, the four-current J, and the derivative 
K of the four-current, as E — E^ip + Ef>K + Ejv 
= Edip + (Wrk + Wjv)/D with Wab defined as 
AB^ — A^B. The term Ejv is zero in case of a 
longitudinally moving point charge, and non-zero in 



case of transverse currents generated by the mag- 
netic field. This contribution therefore vanishes for 
shower trajectories parallel to the geo-magnetic field, 
and the contribution increases with increasing angle 
between shower axis and magnetic field. The most 
important contribution though is the contribution 
Erk caused by the time variation K of the four- 
current. Here we can identify actually two contribu- 
tions: one due to time variation of the charge excess, 
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Figure 26: The x- (full line), y- (dashed) and z-component (dotted) of the electric field of a shower with 6 = 27", 
= 0° , with the full geo-magnetic field, for four observers south, east, north, and west of the point of impact, at 
a distance of 500 m. The upper four plots show the the fields without the dipole contribution, the lower four plots 
represent all contributions. The vector components refer to the Earth frame. 



which is always present, and an second one due to 
the time variation of the transverse current due to 
the geo-magnetic field. This second contribution is 
only present in case of a non-zero angle between the 
shower axis and the magnetic field. For the shower 
we investigated {9 — 27°, ip — 0°), the second contri- 
bution is four times larger than the first. The field is 
thus clearly dominated by the contributions due to 
the transverse current in the geo-magnetic field, but 



the charge excess contribution can not be neglected. 
Also the dipole field i?dip contributes significantly. 

We get in any case a characteristic bipolar signal, 
refiecting first of all the time dependence of the vari- 
ation K of the currents. In case of geo-magnetic con- 
tributions this bipolarity is even more pronounced, 
since here we get an additional contribution Ejv 
whose time signal refiects the time dependence of 
the current J, giving therefore a much broader sig- 
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nal, of opposite sign compared to Ehk- The sum is 
a bipolar signal with two strong pulses. 

We emphasis that this binolarity is in qualitative 
contrast to the findings of [1, [HI, E3| , where 

a microscopic picture has been employed, leading 
to a single pulse. The fundamental difference is that 
the time variation of the currents in the macroscopic 
picture treats in an effective way the creation and 
disappearance of charges. In a microscopic picture 
this has to be treated explicitly, otherwise one misses 
an important part. 



Appendix A: POTENTIAL AND FIELD 
FROM A POINTLIKE CURRENT 



We have 

= / rfa:'°J'3(:r'0),5(4)(a:-^(a;'")), (Al) 



with x° = ct and ^0(2;'°) = x'^ ■ The potential due 
to the current J is 

A^{x) = /io 1 d%^e{x" ^ y'>)Si{x - yf) (A2) 

X J da;'°J^(x"')<5(4)(y-e(a^'°))- 
We change / d'^y... J dx'^ ... into 

dx">... I d^y5\y~i{x"^))..., (A3) 
which gives 

AP[x) = y" (A4) 

x5{{x-ax")f)j'{^"). 
Using R — X — ^{x'^)^ we have 

Af^{x) = ^ / dx"' e{x'' ~ e{x"'))5{R^)jf'{x"'). 
2vr J 

(A5) 
(A6) 



Using 



2, _ 5{x'° - x°) 



{dRydx'^)^,o^^o I ' 



with x*° being the root of R^ (supposing there is only 
one), and 

dR 

d^ = ^^d^ - 
with V = d^/dx"^, we get for xq ^ ^o{x'°) 



Mo J" 



47r \RV\'' 



(A7) 



with all quantities considered at the retarded time 
rjj^g electric field is given as 



E' = c - 



We have 



47r \RV\' 



(A8) 



(A9) 



We have 



_ dJ^/di° 



\RV\ \RV\ " RV\RV 

The derivative of RV is 

d'^{RV] -VV d^x*^ . 

From i?^ = we conclude 

= 2R" - 2RVd"x*° = 0, 

so 



d°'x*° 



RV 



Putting all together, we obtain 



Aw \RV\^ ' 



■d^lRV}. 
(AlO) 

(All) 
(A12) 

(A13) 
(A14) 



with 



L°''^ = RV R^jf^ - RV J'^V"' + VV R°'J'^. (A15) 

where the "dot" means derivative with respect to 
a;*". The electric field is then 



E\x) 



1 L'" - 



Aireoc \RV\^ 



(A16) 



with n being the index of refraction at the observer 
position. 



Appendix B: CURRENTS FROM AIR 
SHOWER SIMULATIONS 

We compute the current by using the CONEX 
shower simulation program, as 

J= ev,~ ^ evj, (Bl) 

positrons! electrons j 

considering all electrons and positrons present at a 
given time. With being the creation time of a 
particle, its velocity some time later can be written 
as 

v{t) — v\\ {to) + |'!7L(to)| (ei cosoji + 6*2 sinw^) , 

(B2) 
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where v\\ and iJj^are the velocity components parallel 
and orthogonal to the magnetic field, uj = ±eB/jm, 
and where {ci} is a system of orthonormal vectors, 
defined via 

e3 = B/\B\, ei ^ vj_{to)/\vj_{to)\, 62 = 6*1X63. 

(B3) 

We define a "shower frame" {A, Ux, Uy, Wz}, where 
the vectors can be expressed with respect to the 
Earth frame basis as 

/ \ / cos 6* \ / sin 6* \ 

= 1 , Uy = , wz - . 
\0 J \ sine J \ - cose J 

(B4) 

The vector u^is parallel to the shower axis (z is the 
longitudinal shower variable). Again with respect 
to the Earth frame basis, the magnetic field may be 
written as 

(sin a cos (f \ 
sin a sin (y5 , (B5) 
cos a / 

where a is the angle between B and the Wz, and (f 
the angle between the shower plane and the plane 
containing B and the vertical axis. Computing 



Bux.y.z, we obtain B in the shower coordinate sys- 
tem, and as well 63 — B/B, 

(sin a sin (p \ 
sin a cos cos 6* + cos a sin 6* . (B6) 
sin a cos (/) sin 6* — cos a cos 9 J 

Then 

v\\ = ■1)63, iT|| = U||e3, , vj_ = v-vi\, v± = \v±\, 

(B7) 

and 



ei^v±/v±, 6*2 = 61X63. (B8) 

Now we know the components of the vectors e*i in 
the shower frame, which allows us to evaluate eq. 
IB2I and determine the components J* of the current 
in the shower frame. We finally obtain easily the 
components of the currents with respect to the Earth 
frame basis as 

J^J'^Ux + jyUy + J'^Uz, (B9) 

and its derivatives. 
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